In geodesy and navigation, a meridian arc is the curve between two points near the Earth's surface having the same longitude. The term may refer either to a segment of the meridian, or to its Arc length. Both the practical determination of meridian arcs (employing measuring instruments in field campaigns) as well as its theoretical calculation (based on geometry and abstract mathematics) have been pursued for many years.
The earliest determinations of the size of a spherical Earth required a single arc. Accurate survey work beginning in the 19th century required several arc measurements in the region the survey was to be conducted, leading to a proliferation of reference ellipsoids around the world. The latest determinations use astro-geodetic measurements and the methods of satellite geodesy to determine reference ellipsoids, especially the geocentric ellipsoids now used for global coordinate systems such as WGS 84 (see numerical expressions).
In 1687, Isaac Newton had published in the Principia as a proof that the Earth was an oblate spheroid of flattening equal to .Isaac Newton: Principia, Book III, Proposition XIX, Problem III, translated into English by Andrew Motte. A searchable modern translation is available at 17centurymaths. Search the following pdf file for 'spheroid'. This was disputed by some, but not all, French scientists. A meridian arc of Jean Picard was extended to a longer arc by Giovanni Domenico Cassini and his son Jacques Cassini over the period 1684–1718.. Freely available online at Archive.org and Forgotten Books (). In addition the book has been reprinted by Nabu Press (), the first chapter covers the history of early surveys. The arc was measured with at least three latitude determinations, so they were able to deduce mean curvatures for the northern and southern halves of the arc, allowing a determination of the overall shape. The results indicated that the Earth was a prolate spheroid (with an equatorial radius less than the polar radius). To resolve the issue, the French Academy of Sciences (1735) undertook expeditions to Peru (Pierre Bouguer, Louis Godin, de La Condamine, Antonio de Ulloa, Jorge Juan) and to Lapland (Maupertuis, Alexis Clairaut, Camus, Le Monnier, Abbe Outhier, Anders Celsius). The resulting measurements at equatorial and polar latitudes confirmed that the Earth was best modelled by an oblate spheroid, supporting Newton. However, by 1743, Clairaut's theorem had completely supplanted Newton's approach.
By the end of the century, Jean Baptiste Joseph Delambre had remeasured and extended the French arc from Dunkirk to the Mediterranean Sea (the meridian arc of Delambre and Méchain). It was divided into five parts by four intermediate determinations of latitude. By combining the measurements together with those for the arc of Peru,
ellipsoid shape parameters were determined and the distance between the Equator and pole along the Paris Meridian was calculated as as specified by the standard toise bar in Paris. Defining this distance as exactly led to the construction of a new standard metre bar as toises.
Apart from the obvious consideration of safe access for French surveyors, the Paris meridian was also a sound choice for scientific reasons: a portion of the quadrant from Dunkirk to Barcelona (about 1000 km, or one-tenth of the total) could be surveyed with start- and end-points at sea level, and that portion was roughly in the middle of the quadrant, where the effects of the Earth's oblateness were expected not to have to be accounted for.
The expedition would take place after the Anglo-French Survey, thus the French meridian arc, which would extend northwards across the United Kingdom, would also extend southwards to Barcelona, later to Balearic Islands. Jean-Baptiste Biot and François Arago would publish in 1821 their observations completing those of Delambre and Mechain. It was an account of the length's variations of portions of one degree of amplitude of the meridian arc along the Paris meridian as well as the account of the variation of the seconds pendulum's length along the same meridian between Shetland and the Balearic Islands.Jean-Jacques Levallois, La méridienne de Dunkerque à Barcelone et la détermination du mètre (1792–1799), Vermessung, Photogrammetrie, Kulturtechnik, 89 (1991), 375–380.
The task of surveying the meridian arc fell to Pierre Méchain and Jean-Baptiste Delambre, and took more than six years (1792–1798). The technical difficulties were not the only problems the surveyors had to face in the convulsed period of the aftermath of the Revolution: Méchain and Delambre, and later François Arago, were imprisoned several times during their surveys, and Méchain died in 1804 of yellow fever, which he contracted while trying to improve his original results in northern Spain.
The project was split into two parts – the northern section of 742.7 km from the belfry of the Church of Saint-Éloi, Dunkirk to Rodez Cathedral which was surveyed by Delambre and the southern section of 333.0 km from Rodez to the Montjuïc Fortress, Barcelona which was surveyed by Méchain. Although Méchain's sector was half the length of Delambre, it included the Pyrenees and hitherto unsurveyed parts of Spain.
Delambre measured a baseline of about 10 km (6,075.90 toises) in length along a straight road between Melun and Lieusaint. In an operation taking six weeks, the baseline was accurately measured using four platinum rods, each of length two toises (a toise being about 1.949 m). Thereafter he used, where possible, the triangulation points used by Nicolas Louis de Lacaille in his 1739–1740 survey of Paris meridian from Dunkirk to Collioure. Méchain's baseline was of a similar length (6,006.25 toises), and also on a straight section of road between Vernet (in the Perpignan area) and Salces (now Salses-le-Château).
To put into practice the decision taken by the National Convention, on 1 August 1793, to disseminate the new units of the decimal metric system, it was decided to establish the length of the metre based on a fraction of the meridian in the process of being measured. The decision was taken to fix the length of a provisional metre (French: mètre provisoire) determined by the measurement of the Paris meridian from Dunkirk to Collioure, which, in 1740, had been carried out by Nicolas Louis de Lacaille and Cesar-François Cassini de Thury. The length of the metre was established, in relation to the toise of the Academy also called toise of Peru, at 3 feet 11.44 lines, taken at 13 degrees of the temperature scale of René-Antoine Ferchault de Réaumur in use at the time. This value was set by legislation on 7 April 1795. It was therefore metal bars of 443.44 that were distributed in France in 1795-1796. This was the metre installed under the arcades of the rue de Vaugirard, almost opposite the entrance to the Senate.
End of November 1798, Delambre and Méchain returned to Paris with their data, having completed the survey to meet a foreign commission composed of representatives of Batavian Republic: Henricus Aeneae and Jean Henri van Swinden, Cisalpine Republic: Lorenzo Mascheroni, Kingdom of Denmark: Thomas Bugge, Kingdom of Spain: Gabriel Císcar and Agustín de Pedrayes, Helvetic Republic: Johann Georg Tralles, Ligurian Republic: Ambrogio Multedo, Kingdom of Sardinia: Prospero Balbo, Antonio Vassali Eandi, Roman Republic: Pietro Franchini, Tuscan Republic: Giovanni Fabbroni who had been invited by Talleyrand. The French commission comprised Jean-Charles de Borda, Barnabé Brisson, Charles-Augustin de Coulomb, Jean Darcet, René Just Haüy, Joseph-Louis Lagrange, Pierre- Simon Laplace, Louis Lefèvre-Ginneau, Pierre Méchain and Gaspar de Prony.
In 1799, a commission including Johann Georg Tralles, Jean Henri van Swinden, Adrien-Marie Legendre, Pierre-Simon Laplace, Gabriel Císcar, Pierre Méchain and Jean-Baptiste Delambre calculated the distance from Dunkirk to Barcelona using the data of the triangulation between these two towns and determined the portion of the distance from the North Pole to the Equator it represented. Pierre Méchain's and Jean-Baptiste Delambre's measurements were combined with the results of the French Geodetic Mission to the Equator and a value of was found for the Earth's flattening. Pierre-Simon Laplace originally hoped to figure out the Earth ellipsoid problem from the sole measurement of the arc from Dunkirk to Barcelona, but this portion of the meridian arc led for the flattening to the value of considered as unacceptable. This value was the result of a conjecture based on too limited data. Another flattening of the Earth was calculated by Delambre, who also excluded the results of the French Geodetic Mission to Lapland and found a value close to combining the results of Delambre and Méchain arc measurement with those of the Spanish-French Geodetic Mission taking in account a correction of the astronomic arc. The distance from the North Pole to the Equator was then extrapolated from the measurement of the Paris meridian arc between Dunkirk and Barcelona and was determined as toises. As the metre had to be equal to one ten-millionth of this distance, it was defined as 0.513074 toise or 3 feet and 11.296 lines of the Toise of Peru, which had been constructed in 1735 for the French Geodesic Mission to Peru. When the final result was known, a bar whose length was closest to the meridional definition of the metre was selected and placed in the National Archives on 22 June 1799 (4 messidor An VII in the Republican calendar) as a permanent record of the result.
For longer arcs, the length follows from the subtraction of two meridian distances, the distance from the equator to a point at a latitude .
This is an important problem in the theory of map projections, particularly the transverse Mercator projection.
The main ellipsoidal parameters are, , , , but in theoretical work it is useful to define extra parameters, particularly the eccentricity, , and the third flattening . Only two of these parameters are independent and there are many relations between them:
The arc length of an infinitesimal element of the meridian is (with in radians). Therefore, the meridian distance from the equator to latitude is
Even though latitude is normally confined to the range , all the formulae given here apply to measuring distance around the complete meridian ellipse (including the anti-meridian). Thus the ranges of , , and the rectifying latitude , are unrestricted.
The calculation (to arbitrary precision) of the elliptic integrals and approximations are also discussed in the NIST handbook. These functions are also implemented in computer algebra programs such as Mathematica Mathematica guide: Elliptic Integrals and Maxima. Maxima, 2009, A computer algebra system, version 5.20.1.
Richard Rapp gives a detailed derivation of this result.Rapp, R, (1991), §3.6, pp. 36–40.
In 1837, Friedrich Bessel obtained one such series, which was put into a simpler form by Helmert,Helmert, F. R. (1880): Die mathematischen und physikalischen Theorieen der höheren Geodäsie, Einleitung und 1 Teil, Druck und Verlag von B. G. Teubner, Leipzig, § 1.7, pp. 44–48. English translation (by the Aeronautical Chart and Information Center, St. Louis) available at Krüger, L. (1912): Konforme Abbildung des Erdellipsoids in der Ebene. Royal Prussian Geodetic Institute, New Series 52, page 12
Because changes sign when and are interchanged, and because the initial factor is constant under this interchange, half the terms in the expansions of vanish.
The series can be expressed with either or as the initial factor by writing, for example,
Because this series provides an expansion for the elliptic integral of the second kind, it can be used to write the arc length in terms of the geodetic latitude as
Delambre and Bessel both wrote their series in a form that allows them to be generalized to arbitrary order. The coefficients in Bessel's series can expressed particularly simply
and is the double factorial, extended to negative values via the recursion relation: and .
The coefficients in Helmert's series can similarly be expressed generally by
This result was conjectured by Friedrich HelmertHelmert (1880), §1.11 and proved by Kazushige Kawase.Kawase, K. (2011): A General Formula for Calculating Meridian Arc Length and its Application to Coordinate Conversion in the Gauss-Krüger Projection, Bulletin of the Geospatial Information Authority of Japan, 59, 1–13
The extra factor originates from the additional expansion of appearing in the above formula and results in poorer convergence of the series in terms of compared to the one in .
Substituting the values for the semi-major axis and eccentricity of the WGS84 ellipsoid gives
On the ellipsoid the exact distance between parallels at and is . For WGS84 an approximate expression for the distance between the two parallels at ±0.5° from the circle at latitude is given by
The quarter meridian can be expressed in terms of the complete elliptic integral of the second kind,
The quarter meridian is also given by the following generalized series:
The numerical expression for the quarter meridian on the WGS84 ellipsoid is
Alternatively, Helmert's series for the meridian distance can be reverted to giveHelmert (1880), §1.10Adams, Oscar S (1921). Latitude Developments Connected With Geodesy and Cartography. US Coast and Geodetic Survey Special Publication No. 67. p. 127.
Similarly, Bessel's series for in terms of can be reverted to giveHelmert (1880), §5.6
Adrien-Marie Legendre showed that the distance along a geodesic on a spheroid is the same as the distance along the perimeter of an ellipse. For this reason, the expression for in terms of and its inverse given above play a key role in the solution of the geodesic problem with replaced by , the distance along the geodesic, and replaced by , the arc length on the auxiliary sphere.Helmert (1880), Chap. 5 The requisite series extended to sixth order are given by Charles Karney, Addenda. Eqs. (17) & (21), with playing the role of and playing the role of .
19th century
19th century
The nautical mile
Calculation
f&=\frac{a-b}{a}\,, \qquad e^2=f(2-f)\,, \qquad n=\frac{a-b}{a+b}=\frac{f}{2-f}\,,\\
b&=a(1-f)=a\sqrt{1-e^2}\,,\qquad e^2=\frac{4n}{(1+n)^2}\,.
\end{align}
Definition
m(\varphi) &=\int_0^\varphi M(\varphi) \, d\varphi \\
&= a(1 - e^2)\int_0^\varphi \left(1 - e^2 \sin^2 \varphi \right)^{-\frac32} \, d\varphi\,.
\end{align}
The distance formula is simpler when written in terms of the
parametric latitude,
where and .
Relation to elliptic integrals
It may also be written in terms of incomplete elliptic integrals of the second kind (See the NIST handbook Section 19.6(iv)),
m(\varphi) &= a\left(E(\varphi,e)-\frac{e^2\sin\varphi\cos\varphi}{\sqrt{1-e^2\sin{}^{\!2}\varphi}}\right) \\
&= a\left(E(\varphi,e)+\frac{d^2}{d\varphi^2}E(\varphi,e)\right) \\
&= b E(\beta, ie')\,.
\end{align}
Series expansions
Expansions in the eccentricity ()
where
D_0 &= 1 + \tfrac{3}{4} e^2 + \tfrac{45}{64} e^4 + \tfrac{175}{256} e^6 + \tfrac{11025}{16384} e^8 + \cdots, \\5mu
D_2 &= - \tfrac{3}{8} e^2 - \tfrac{15}{32} e^4 - \tfrac{525}{1024} e^6 - \tfrac{2205}{4096} e^8 - \cdots, \\5mu
D_4 &= \tfrac{15}{256} e^4 + \tfrac{105}{1024} e^6 + \tfrac{2205}{16384} e^8 + \cdots, \\5mu
D_6 &= - \tfrac{35}{3072} e^6 - \tfrac{105}{4096} e^8 - \cdots, \\5mu
D_8 &= \tfrac{315}{131072} e^8 + \cdots.
\end{align}
Expansions in the third flattening ()
with
H_0 &= 1 + \tfrac{1}{4} n^2 + \tfrac{1}{64} n^4 + \cdots, \\
H_2 &= - \tfrac{3}{2} n + \tfrac{3}{16} n^3 + \cdots,&
H_6 &= - \tfrac{35}{48} n^3 + \cdots, \\
H_4 &= \tfrac{15}{16} n^2 - \tfrac{15}{64} n^4 - \cdots,\qquad&
H_8 &= \tfrac{315}{512} n^4 - \cdots.
\end{align}
and expanding the result as a series in . Even though this results in more slowly converging series, such series are used in the specification for the transverse Mercator projection by the National Geospatial-Intelligence AgencyJ. W. Hager, J.F. Behensky, and B.W. Drew, 1989. Defense Mapping Agency Technical Report TM 8358.2. The universal grids: Universal Transverse Mercator (UTM) and Universal Polar Stereographic (UPS) and the Ordnance Survey of Great Britain. A guide to coordinate systems in Great Britain, Ordnance Survey of Great Britain.
Series in terms of the parametric latitude
m(\varphi) = \frac{a+b}2\left(B_0\beta + B_2\sin2\beta + B_4\sin4\beta + B_6\sin6\beta + B_8\sin8\beta + \cdots\right),
with
B_0 &= 1 + \tfrac{1}{4} n^2 + \tfrac{1}{64} n^4 + \cdots = H_0\,,\\
B_2 &= - \tfrac{1}{2} n + \tfrac{1}{16} n^3 + \cdots, &
B_6 &= - \tfrac{1}{48} n^3 + \cdots, \\
B_4 &= - \tfrac{1}{16} n^2 + \tfrac{1}{64} n^4 + \cdots, \qquad&
B_8 &= - \tfrac{5}{512} n^4 + \cdots.
\end{align}
m(\varphi) = \frac{a+b}2\Biggl(
&B_0\varphi - B_2\sin2\varphi + B_4\sin4\varphi - B_6\sin6\varphi + B_8\sin8\varphi-\cdots \\[-3mu]
&\quad -\frac{2n\sin2\varphi}{\sqrt{1+2n\cos2\varphi+n^2}}\Biggr).
\end{align}
Generalized series
\begin{cases} c_0\,, & \text{if }k = 0\,, \\5px
\dfrac{c_k}{k}\,, & \text{if } k > 0\,,
\end{cases}
where
Numerical expressions
m(\varphi)&=\left(111\,132.952\,55\,\varphi^{(\circ)}-16\,038.509\,\sin 2\varphi+16.833\,\sin4\varphi-0.022\,\sin6\varphi+0.000\,03\,\sin8\varphi\right)\mbox{ metres} \\
&= \left(111\,132.952\,55\,\beta^{(\circ)}-5\,346.170\,\sin 2\beta-1.122\,\sin4\beta-0.001\,\sin6\beta-0.5\times10^{-6}\,\sin8\beta\right)\mbox{ metres,}
\end{align}
where is expressed in degrees (and similarly for ).
Quarter meridian
It was part of the historical definition of the metre and of the nautical mile, and used in the definition of the hebdomometre.
where are the first and second eccentricities.
(For the formula of c0, see section #Generalized series above.)
This result was first obtained by James Ivory.
m_\mathrm{p} &= 0.9983242984312529\ \frac{\pi}{2}\ a\\
&= 10\,001\,965.729\mbox{ m.}
\end{align}
Full meridian (polar perimeter)
The perimeter of a meridian ellipse can also be rewritten in the form of a rectifying circle perimeter, . Therefore, the rectifying Earth radius is:
It can be evaluated as .
The inverse meridian problem for the ellipsoid
until convergence. A suitable starting guess is given by where
is the rectifying latitude. Note that it there is no need to differentiate the series for , since the formula for the meridian radius of curvature can be used instead.
where
H'_2 &= \tfrac{3}{2} n - \tfrac{27}{32} n^3 + \cdots,&
H'_6 &= \tfrac{151}{96} n^3 + \cdots, \\
H'_4 &= \tfrac{21}{16} n^2 - \tfrac{55}{32} n^4 + \cdots,\qquad&
H'_8 &= \tfrac{1097}{512} n^4 + \cdots.
\end{align}
where
B'_2 &= \tfrac{1}{2} n - \tfrac{9}{32} n^3 + \cdots,&
B'_6 &= \tfrac{29}{96} n^3 - \cdots, \\
B'_4 &= \tfrac{5}{16} n^2 - \tfrac{37}{96} n^4 + \cdots,\qquad&
B'_8 &= \tfrac{539}{1536} n^4 - \cdots.
\end{align}
See also
External links
|
|